Effective electron-electron interaction in a two dimensional paramagnetic system 
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We analyze the effective electron-electron interaction in a two dimensional polarized paramagnetic 
system. The spin degree of freedom, s, is manifestly present in the expressions of spin dependent 
local field factors that describe the short range exchange (x) and correlation (c) effects. Starting from 
the exact asymptotic values of the local field correction functions for large and small momentum at 
zero frequency we obtain self-consistent expressions across the whole spectrum of momenta. Then, 
the effective interaction between two electrons with spins s and s' is calculated. We find that the 
four effective interactions, up-up, up-down, down-up and down-down, are different. We also obtain 
their qualitative dependence on the electronic density and polarization and note that these results 
are independent of the approximation used for the local field correction functions at intermediate 
momenta. 

PACS numbers: 71.10.-w,72.25.-b,71.45.Gm 
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I. INTRODUCTION 

To accomplish spin dependent conduction in electronic, 
devices has become a very intense quest in recent years.u 
The II- VI dilute magnetic semiconductors, like CdMnTc 
and ZnMnSe, seem to be the most promising materials 
for achieving this goal. At low doping levels these sys- 
tems exhibit an enhanced Zeeman splitting of the elec- 
tronic levels that arises from strong Kondo-like exchange 
between the embedded Mn ions and the delocalized elec- 
tronic states. Moreover, the II- VI dilute magnetic semi- 
conductors are n-type conductors with relative good mo- 
bility and large spin coherence times (~ 10ns). However, 
an open question remains as to how the magnetic interac- 
tion between the localized spins and the itinerant carriers 
reflects on the transport properties of these structures. 
The purpose of this work is to analyze some aspects of 
this problem. 

Since the Zeeman splitting of the electron levels 
was found jtp be independent of the local magnetic 
environments, a simple model of a II- VI magnetic semi- 
conductor is a spin polarized electron gas in the presence 
of a static magnetic field. The field lifts the spin de- 
generacy and induces an equilibrium polarization, £ = 
(n-f — n^)/(rif +JijJ. The strong magnetic interaction be- 
tween the itinerant carriers and the localized spins is re- 
flected in the large value of the effective gyromagnetic fac- 
tor, 7*, up to hundreds of times the band value. On ac- 
count of the large 7* , even low magnetic fields are enough 
to produce large polarizations, and it is assumed that £ 
can vary continuously between — 1 and 1 as a function of 
the static magnetic field. This approximation integrates 
out the degrees of freedom of the static spins under the 
enlarged value of 7* and focuses on the itinerant carriers 
and the many body interaction among them. The lat- 
ter is independent of the source of the spin polarization, 
and we expect our results to maintain their validity also 
in the case of a self-consistent magnetic field as source 
of the spin polarization, a situation which is consistent 
with an itinerant ferromagnet. The model can also be 



extended to describe the paramagnetic state of the III-V 
magnetic heterostructures where a mostly uniform inter- 
nal magnetic field is created by having a magnetic 
density much larger than the itinerant-carrier density. 
However, the small electronic densities in the III-V based 
compounds make difficult the blind use of the paramag- 
netic electron gas model. 

The explicit spin dependence of the electron-electron 
interaction becomes manifest when the exchange (x) and 
correlation (c) effects of the local Coulomb repulsion are 
included. At finite values of the polarization, the many- 
body short range interaction, which is density dependent, 
is different for up and down spin electrons. A realis- 
tic picture of the electronic interaction and its screen- 
ing is obtained by using spin dependent local field cor- 
rection functions, G^' c (q, oj), that describe the exchange 
and correlation hole around each electron. The self con- 
sistent treatment of exchange and correlation effects has 
proved very important in understanding the physics of 
normal metalsjj but to our knowledge it has not been 
fully analyzed in spin-polarized systems, where the spin 
dependence of the local field corrections becomes man- 
ifest. In addition, the relevance of exchange and corre- 
lation effects increases as the dimensionality of the elec- 
tron gas is lowered. For example, the importance of lo- 
cal effects is clearly reflected in the dielectric function 
of an unpolarized two-dimensional electron gas, which 
becomes overscreened in a wide range of momenta and 
electronic density!] However, calculations where these ef- 
fects are absent, such as the conventional random phase 
approximation, predict a positive dielectric function for 
the whole parameter range. Therefore, we believe it is 
fundamental to include these local effects in a complete 
treatment of the quasi-two-dimensional diluted magnetic 
semiconductors. 

Obtaining the exact frequency and wave vector de- 
pendence of the local field corrections is a very difficult 
problem which remains unsolved even in the case of the 
unpolarized electron system. Fortunately, the asymp- 
totic values of the local field factors can be obtained 
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exactly in two limiting cases. At zero frequency and 
small wavevectors, sum rules are used to connect the 
static limits of the response functions to certain ther- 
modynamic coefficients. For large frequency and large 
wavevector, an iterative method generates the exact e: 
pressions for the local field functions up to second order 
This approach uses the equation of motion satisfied by 
the Wigner distribution function of the particle density. 

Numerical estimates of the response functions of the 
three dimensional unpolarized electron gas have shown 
that local field factors smoothly interpolate between the 
asymptotic small and large wave-vector behaviorB This 
feature is expected to exist also in the case of a spin polar- 
ized system, and, consequently, we use the exact asymp- 
totic values of G^Lclmj) for large and small momentum 
at zero frequency0E3E!J as a starting point in obtaining 
their approximate expressions across the whole spectrum 
of momentum. Then, G a (q, lo) are used to calculate the 
spin dependent effective electron interaction potentials, 
Hrer' (q, w). We compare two interpolation schemes to 
conclude that the effective potentials are qualitatively in- 
dependent of the approximation used for the local field 
correction functions at intermediate momenta. 

In section II, we present the self consistent formula- 
tion of the effective interaction that incorporates the spin 
dependent local field correctiaas-ibllowing the approach 
of Kukkonen and Overhausprc3ll3 generalized for a spin 
polarized electron systemO This formalism permits the 
derivation of the effective spin dependent scattering po- 
tentials experienced by an electron and of the response 
functions. In section III, we discuss in detail the proce- 
dure used to obtain the approximate expressions for the 
local field factors. Section IV presents our conclusions. 
The appendix is dedicated to a derivation of the corre- 
lation function of a two dimensional electron gas at the 
origin. 



II. EFFECTIVE INTERACTION POTENTIALS 

The simplest approximation of the effective interac- 
tion in the electron gas is the random phase approxi- 
mation (RPA), which incorporates the screening but ig- 
nores the exchange and correlation effects. Therefore, 
all spin dependence is lost. A more realistic description, 
which considers the short range Coulomb interaction ef- 
fects, was proposed by Hubbardjl3 whose expression for 
the dielectric constant introduces the local field correc- 
tion, a wavevector dependent function that describes the 
difference between the real particle density and its mean 
field (RPA) counterpart. 

The microscopic origin of the local field corrections 
was elucidated by Kukkonen and OverhauserJ13 In their 
model of the electron-electron interaction, the exchange 
(x) and correlation (c) are explicitly incorporated by con- 
sidering associated local field factors, G(q, ui), in the ex- 
pression of the effective potential experienced by one elec- 
tron in the presence of all the rest. In a self-consistent 



V, rat. fit 
V, exp. fit 
J , rat. fit 
J , exp. fit 




FIG. 1: Normalized static effective interactions V{(\,lu — 0) 
and J(q, u = 0) for an unpolarized electron gas (£ = 0) as 
functions of the normalized momentum q/fcf for r s = 4. Re- 
sults using a rational (Eq. (p"l|)) and an exponential (Eq. (p^)) 
fit for the local field correction functions are displayed. Note 
the similarity of the results for the two fittings. Results for 
the random phase approximation are also displayed. 
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FIG. 2: Normalized static effective interactions V" CT (q,aj = 0) 
and J CT (q, u = 0) for a polarization of C = 0.5 as functions 
of the normalized momentum q/kp for r s — 4. Results for 
the random phase approximation are also displayed. We have 
used a rational fit (Eq. (|l l\j ) for the local field corrections. 
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FIG. 3: Normalized effective interaction V e //(q) for an un- 
polarized electron gas (£ = 0) as a function of the normal- 
ized momentum q/fcj- for r s = 4. Results using a rational 
(Eq. ([ll])) and an exponential (Eq. (|l^)) fit for the local field 
correction functions as well as results for the random phase 
approximation are displayed. Note that the RPA interaction 
is the same as that in Fig. |l| 
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FIG. 4: Normalized effective interaction V e //(q) for a polar- 
ization of £ = 0.5 as a function of the normalized momentum 
q/kF for r s = 4. Results using a rational (Eq. (0)) and an 
exponential (Eq. ([[2])) fit for the local field correction func- 
tions as well as results for the random phase approximation 
are displayed. The RPA curve is the same as that in Fig. ^| 



formulation, the details of which will be given below in 
the case of a spin polarized electron system, this theory 
leads to the Hubbard dielectric function. The method is 
equivalent to considering the many body interaction in 
all orders, as it has been shown by using diagrammatic 
techniques .E3 

In the case of a spin polarized system, it is appropri- 
ate to introduce spin dependent local field corrections 
because the many body effects are density dependent, 
and they depend on the spin directly through the exclu- 
sion principle. In a direct generalization of Ref. [l2| the 
effective potential V aa > experienced by an electron with 
spin a in the presence of an second electron_with spin a' 
and charge density p(q) can be written as:E3 



, (q, w) = «(q) { [1 - G+ (q, w)] (p(q) + An(q, w)) 

-otf-(q,w)(tr / p(q) + Am,(q,w))} (1) 

where An = An-p + Ani and Am z = An-t — Ani are the 
particle and spin density fluctuations, respectively, in- 
duced by the presence of the second electron, and v(q) is 
the bare Coulomb interaction, which is equal to 2we 2 /q 
in two dimensions. The local field functions G^(q, uS) 
incorporate the exchange {G%) and correlation effects 
(G% al G% s ), which induce a local decrease in the elec- 
tronic density around an electron of given spin a com- 
pared with its RPA value. G>(q, ui) is the sum of all 
same and opposite spin effects, and G~ (q, to) is the dif- 
ference of same and opposite spin effects: Gr*(q,w) = 
G x a + G c aa ± G c aS JB In brief, Eq. (@) shows how the elec- 
tron charge density is decreased from the corresponding 
RPA value to account for its own short range effects. 

The role played by G+ and G~ in the physical proper- 
ties of the electron gas is quite different. For an unpolar- 
ized electron gas, G + enters the electrical response func- 
tions, and G~ the magnetic response. It is well known 
that in this limit the dielectric function and the electronic 
susceptibility can be written down as functions of G + (q) 
and G~ (q) respectively:E3 



e(q,w) = 1 - 

x(q,^) = 



u 9 II(q,cj) 



l + G+(q)D,n(q, W ) 

- 7 * 2 n(q,^) 
l + G-(q)«,n(q,w) 



(2) 
(3) 



where 7* is the enhanced effective gyromagnetic ratio and 
II = n-j-f + Tin i s the polarization function. However, in 
a polarized gas both local field factors are interconnected 
and appear in the expressions of the dielectric function 
and the magnetic susceptibility. 

In a linear approximation, the density fluctuations are 
proportional to the effective potentials, where the pro- 
portionality coefficients are appropriately defined polar- 
ization functions: An-f = Hfi-Vf CT ' (An^ = n^V^o-'). 
Here II ffff / is the generalized polarization bubble for the 
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non-interacting electron gas, 

iW(q,w) = - 7^ — ( 4 ) 

K k + ^ kcr _ Sk+qo-' 

where £k<r = e k + sign(tr)7*i? is the single particle energy 
in the static magnetic field B, rik, ff is the occupation func- 
tion, and V is the volume of the system. Eq. (Q) gener- 
ates the well known expressions for the polarization func- 
tion of the non-interacting electron gas when the usual 
simplifications are considered: zero temperature and a 
parabolic energy dispersion, ek = fi 2 fc 2 /2m*, where m* 
is the electronic band mass. Neither approximation is 
expected to modify our results. 

Therefore, in two dimensions thp polarization operator 
for the complex frequency iuj isj!3 



K 2qk Fa - 



where q = q/kp is the normalized momentum, kp a = 
kFrjjkp the normalized Fermi momenta of the spin a 



electronic population, and uj = ui/Ep the normalized 
frequency. The use of normalized variables allows us to 
express easily all our results in terms of our only free 
parameter: the effective electronic density, or the ratio 
between r s and the-.effective Bohr radius of the system, 
a* B = h 2 /(m*e* 2 ).Bl The retarded polarizability can be 
obtained by analytical continuation. 

The truly effective interaction potentials, which are 
used in the calculation of scattering amplitudes, are ob- 
tained from V aa ' by subtracting the direct exchange 
and correlation between the two electrons: the term 
v(q)_p(q)[-G+(q) - aa'G-(q)} in Eq. | Thus, the ef- 
fective interaction can be expressed as: 



V<r,a< = V a > - J a >(T ■ a' 



(6) 



where the first term includes the bare Coulomb interac- 
tion and the interaction mediated by charge fluctuations, 
and the second term reflects the interaction mediated by 
spin fluctuations. In an unpolarized electron gas, V and 
J are spin independent. In addition, V depends only 
on G + while J depends only on G~, as it was shown in 
Ref. [l2| However, the imbalance between the two spin 
populations in a spin polarized electron gas induces a 
dependence on the spin index. Thus, V a and J a are: 



V a = [v< l p Ki /D(q,uj)}{l + -[G+ + G± + G- - G^Dfau) -l]+v^[G-(l- G%) + G=(1-G+)]} (7) 
J, = -KpJD{< l ,u)][±[G+ -Gt + G- +Gr][ J D(q,o;)-l]+ Uq n^[G-(l-G±) + Gr(l-G+)]} (8) 
where £>(q,w) = [1 - « q (l - G+ - Gf )n TT ][l - « q (l - G + - G7/)II U ] - [1 - G+ + G£ ][1 - G+ + G x > q n TT n u . 



Figs. |lj and g show the static effective interactions 
V^(q,w 0) and J CT (q, u> = 0) between point-like 
electrons^! in a two dimensional electron gas as func- 
tions of the normalized momentum q/kp. Fig. |l] displays 
the results for an unpolarized electron gas. Fig. || corre- 
sponds to a £ = 0.5 polarization. We have chosen r s = 4, 
where r s is measured in units of the effective Bohr ra- 
dius of the systemO The expressions used for the local 
field corrections are discussed in the following section. 
For comparison, we also include the RPA effective inter- 
action potential, VppA(q,uj — 0), obtained from equa- 
tion (0) by neglecting the local field corrections. Note 
that Jrpa(<1, uj = 0) = 0, since the interaction mediated 
by spin fluctuations is directly proportional to the local 
corrections. 

Fig. [l] shows clearly the importance of the local effect 
in the calculation of physical properties. First, as we 
have already mentioned, the spin mediated interaction, 
J<r(q, uj), becomes noticeable. Second, V a (q,uj) is greatly 
enhanced at small momenta in comparison with the RPA 
prediction. In addition, Fig. shows how both effec- 



tive interactions split when the electron gas is polarized, 
an effect also absent in the RPA approach. The split- 
ting of J ff (q, u/) is more pronounced that the splitting of 
V a (q, uj) and, even, becomes negative for intermediate 
momenta. 

Larger values of the polarization induce larger splitting 
of the up and down effective interactions until a spin 
wave instability is reached. For the density value we are 
considering (r s = 4), this instability happens at C ~ 0.72. 
At that value of the polarization, J CT (q, u>) and V a (q,uj) 
diverge for a certain value of the momenta that define 
the periodicity of the spin density wave. 

In order to have a complete picture of the effective in- 
teractions, we also discuss the effective interaction which 
is used to calculate the electronic selfenergy. It can be de- 
rived in a similar manner using the self-consistent relation 
established between the density fluctuations ajad the ef- 
fective potential induces by an external charge.E3 This ef- 
fective interaction can be written as V.ff = — ; — — — r , 

ett e(q,w = 0)' 

where e(q, uj) is the dielectric function of a polarized elec- 



5 



tron gas: 



e(q,w) = l-u q - 



n 
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+ n u +2« q (Gf + Gpn TT n u 



(9) 



where D(q,w) = [1 + v q (G+ + Gf )n TT ][l + w q (G+ 



Gpn u ] - [G+ - Gf][G+ - G^] Uq n TT n u . 



Fig. || displays 14// (q) for an unpolarized electron gas 
with the same density as before (r s = 4). Since V e //(q) 
represents the electrostatic potential seen by a spinless 
test charge J13 it is very illustrative to compare it with 

the RPA result, V RPA = ^ r-@ We find that 

ei?,PA(q, u = 0) 
by including the local factors the effective interaction be- 
comes overscreened for q less than a critical value, q c , 
which depends on the electronic density and the polariza- 
tion. Similar results for zero polarization were reported 
previouslyo The critical value of the momentum q c de- 
pends weakly on the approximation used for the local 
field factors. Therefore, the effective interaction with lo- 
cal field corrections is very different from the RPA effec- 
tive interaction which is always positive. Fig. || displays 
Ve//(q) for an electron gas with £ = 0.5. The over- 
screened region is also present, and the singularities at 
the Fermi momenta of the majority and minority spin 
populations are clearly seen. 
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FIG. 5: Initial slope of the local field correction function otj 
as a function of the polarization (Q for different values of r s 
as indicated in the legend. 



III. SPIN DEPENDENT LOCAL FIELD 
CORRECTIONS 

Any quantitative calculation of the effects of the 
electron-electron interaction on many physical proper- 
ties of interest requires the knowledge of the correct 
form of the local field correction functions. As we have 
already mentioned, obtaining the exact frequency and 
wave vector dependence of these functions has been an 
elusive problem. Our first approximation is to neglect 
the frequency dependence of the local field corrections. 
Although the local field functions represent a dynami- 
cal effect, -they vary slowly on the scale of the Fermi 
frequency,trL3 and it is acceptable to neglect their fre- 
quency dependence. _ . 

The exact asymptotic values of G^(q, lo) for largeO 
and small momentum!] at zero frequency have been ob- 
tained previously and are summarized in Table |, where 
q = q/kp and g(0) is the two-particle correlation function 
at the origin. Note that the values of G*' (q = 0) in- 
cluded in the table have been derived by considering only 
the contribution coming from the electronic exchange en- 
ergy. The proper expression for G^(q = 0) incorporates 
also a combinatioruof derivatives of the correlation en- 
ergy of the system,Q that need to be added to the values 
given in Table [j| 

If only the exchange contribution is considered, 
G + (q = 0) = G~(q = 0) for the unpolarized electron 
gas. Using the asymptotic values from Table 1 and the 



previously obtained expressions for the dielectric func- 
tion [Eq. (||)] and magnetic susceptibility [Eq. (B)], we 
find that in the static limit both response functions de- 
velop a pole for r s > ~ 2.221. Note that a pole 
in the dielectric function does not indicate an instability 
of the paramagnetic phase, it just implies the appear- 
ance of a region in momentum space where the Coulomb 
interaction is overscreenedEj However, a pole in the sus- 
ceptibility points towards a spin density wave instability. 
Since numerical calculations have shqwn the stability of 
the paramagnetic phase for r s < 30E3 it is crucial to in- 
clude the correlation energy in the calculation of the local 
field factors. Including the correlations in G^r prevents 
the occurrence of a spin/charge density wave instability 
in the two-dimensional unpolarized electron gas. 

Even though there are numerous calculations of the 
correlation energy of an unpolarized electron gas, approx- 
imate expressions for the correlation energy of a polarized 
electron gas are less numerous due to the fact that the 
magnetic response functions are computationally harder 
to obtain. lj For the two dimensional electron gas we us© 
the Monte Carlo calculation by Tanatar and Ceperley,Ell 
where the correlation energy is expressed in the form of 
an analytic interpolation formula: E c (r s , = E c (r Sl C, = 
0) + ( 2 (E c (r s ,(= l)-E c (r B ,C = 0)), where E c (r s , ( = 0) 
and E c (r s , £ = 1) are approximated using a Pade scheme. 

By including the contribution from the correlation en- 
ergy, we derived analytical expressions for the initial 
slope of the local field correction functions <x^ = G^ (q — > 
0)/(q) d_1 , where d is the dimension of the system. These 
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FIG. 6: Initial slope of the local field correction function a7 
as a function of the polarization (£) for different values of r s . 

TABLE I: Exact asymptotic values of the local-field factors 
(q, w = 0) of a two dimensional electron gas for large and 
small a at zero frequency. Gt (q — > oo) was calculated in 
Ref. [llf The values of GJ ± (q = 0) are derived from Ref. § 
considering only the contribution coming from the electronic 
exchange energy. 



G^(q — » oo) 


1 - (1 - C).9(0) 


G^(q — > oo) 


(i-C)s(o) 


G*<+(q^0) 


(g/(27T)){CVl + C+(2-C)Vl-C} 


Gf(q-0) 


(g/(27r)){ [(2 + o/vi + CJ + tc/vi - CJ} 



initial slopes are function of the density and the polar- 
ization of the electron gas. The behavior of a + and a~ 
is quite different. Fig. shows a* for a two dimensional 
gas as function of the polarization for several values of r s . 
Note that at is always positive, as it should be, since lo- 
cal effects always decrease the bare electron charge den- 
sity. On the other hand, a^ depends strongly on the 
polarization for any electronic density as it can be seen 
in Fig. [| where a^ is displayed. In our approximation 
a J diverges for a fully polarized electron gas and, as a 
consequence, G^~ (q) becomes a constant and equal to its 
value at q — > oo. In addition, a^ becomes negative in 
diluted systems and small polarizations H This behav- 
ior favors the stability of the paramagnetic phase. Note 
that similar arguments apply to a^ due to the fact that 

af(C)=af(-C). 

General expressions for G^(q) can be interpolated be- 
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FIG. 7: Local field corrections G^(q,cu = 0) versus normal- 
ized momentum q/fcf for a two dimensional unpolarized elec- 
tron gas with r s = 4. Results for the rational fit (Eq.( [ll|)) 
(solid curves) and for the exponential fit (Eq. ( |l^)) (dashed 
curves) are displayed. The inset is a blow-up of the lower 
corner with the G~ functions. 

tween the large and small momentum limits. The sim- 
plest way to interpolate G"t(fl) for a regular metal was 
first suggested by Hubbard :E£l 

G + (q) = r— (10) 

1+ (a/G+(oo)Jq 2 

In a later work, Singwi and collaborators00 remarked 
that their numerical results for G + (q) in an unpolarized 
three dimensional gas could be adequately fitted by the 

simpler expression G + (q) = Afl — e~ Bq ^j, where the 

two fitting parameters are smoothly varying functions of 
the electronic density. They noted that this expression 
fits their results well at small and intermediate values of 
momentum but not at larger values. However, G + (q) 
is relatively unimportant at large q. Similar conclusions 
have been reached for the unpolarized two dimensional 
electron gas,E-3 where the argument of the exponential 
is linear rather than quadratic. However, the possibility 
that G + (q) have a peak near q — 2kpfJ as a residue of 
the sharp peak in the exchange potential, is a long stand- 
ing issue. Recent calculations show that the inclusion of 
short-range correlations has the effect of washing out this 
peakfl 

Approximate expressions for G _ (q) are less numer- 
ous, partly due to the fact that G~ is related to the 



7 




ized momentum q/fcf for a two dimensional electron gas with 
r s — 4 and polarization ( — 0.5. Results for the rational fit 



(Eq.( 0)) 



Gy (dotted curves) and Gf (solid curves) and for 



the exponential fit (Eq. ( |l^)), Gf (dashed) and Gf (long- 
dashed) are displayed. The inset is also a blow-up of the lower 
corner with the G~ functions. 



magnetic response functions, which are computationally 
challengmgES In addition, it is known that for a three 
dimensional gas, G _ (q) can not be a monotonic func- 
tion of momentum because its slope at q = is positive 
but asymptotically approaches a negative value at large 
momentum. Zhu and Overhauserllia suggested a simple 
rational function with a maximum at q = 2kp. 

Our approach shows that the initial slope of the lo- 
cal field correction functions G+(q) is always positive. 
In contrast, the initial slope of G~(q) could be negative 
for diluted systems and small values of the polarization. 
In the large momentum limit, Gj(q — > oo) is always 
positive in two dimensions. However, in three dimen- 
sions Gt(q — > oo) might be negative for £ < —2/3, and 

G7(q — + oo) is negative for all values of the polarization 
if the electron density verifies r s > 1.18.123 

Given the diversity of behaviors in the limiting values 
of the local field correction functions, we consider a gen- 
eral interpolating scheme for all the cases. For a value 
of the polarization £ and of the electron density r s , we 
calculate the product a^G^(q — * oo) where has been 
defined previously. When a^G^(q — > oo) > 0, we use 



the following fitting expressions for the local field factors: 

— , rational fit; (11) 



1+ a?/G?(oo) q*~i 



G±(q) = G±(oo)(l-exp- 



(a±/G± (oc))q d 



exponential fit. (12) 



In the opposite case, a^G*(q — > oo) < we fit the 
local field factor to the simplest rational or exponential 
function with a maximum at q = 2k f'- 



G±(q) = 



I f ( 7 ±/G±(oo))q d + 1 ' 
G£(q) = G£(oo)(l-exp> 



rational fit; (13) 



where 7^ 



a±(d -l)/4 



2 d (a±/G±(oc))-(d+l) 



exponential fit. (14) 



for the rational 



FIG. 8: Local field corrections G^ (q, u) = 0) versus normal- an< ^ To- 



for the exponential fit. 



2d G±(oo) 

Fig. Q shows the momentum dependence of the local 
field correction functions, G^(q, uj = 0), for an unpolar- 
ized two dimensional electron gas. We have chosen an 
intermediate value for the electron density, r s = 4, which 
corresponds to a density of 1.86 -10 10 cm -2 in GaAs. It is 
clear that both rational and exponential fittings show the 
same qualitative behavior, although the exponential fit- 
ting function approaches the large momentum limit faster 
than the rational one. Fig. || displays Gf (q, uj — 0) ver- 
sus normalized momentum for a polarization of C, = 0.5 
and the same electronic density. It is clear how the local 
field factors split in a polarized gas, although their func- 
tional behavior is the same as that for zero polarization. 

Finally, let us mention that the needed expression of 
the two dimensional two-particle correlation function at 
the origin, which is derived in the appendix, is: 



9(0) = 



1 



2(1 + 0.5857r s ) 2 



(15) 



IV. CONCLUSIONS 

Motivated by recent experimental developments in the 
physics of diluted magnetic semiconductors, we devised a 
framework that explicitly incorporates the spin degree of 
freedom of the itinerant carriers. In our model, the static 
spins are integrated out under the enlarged gyromagnetic 
ratio, and the itinerant electrons are treated as a spin 
polarized gas in the presence of a static magnetic field. 
Therefore, we focused on the itinerant carriers and the 
many body interactions among them, which are modeled 
by using spin dependent local field correction factors. 
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We found approximate expressions for the local field 
correction functions in a two dimensional spin polar- 
ized gas by interpolating their exact asymptotic values 
at small and large wavevectors. Our results indicate 
that the overall behavior of the local field corrections, 
(q) , does not depend on the exact interpolation func- 
tion used. Given the density dependence of the local field 
functions, the values derived in the paramagnetic electron 
gas approximation should represent a realistic estimate 
irrespective of the source of the polarizing magnetic field. 

Using the expressions for the local field factors we cal- 
culated the effective potentials between the two spin pop- 
ulations. In contrast with the RPA approach, the inclu- 
sion of the local factors produces a noticeable value of 
the interaction mediated by spin fluctuations. The in- 
teraction mediated by charge fluctuations is also greatly 
enhanced at small momenta. Charge and spin mediated 
interactions split when the electron gas is polarized, and, 
as a consequence, the effective interaction between two 
electrons with spins a and a depends on the value of a 
and a . 

We also calculated the effective screened interaction, 

V e ff = —. — — — r- , for the two approximate expressions 
e(q,u; = 0) 

of the local field correction factors we used. We found 
that both approximations produce qualitatively equiva- 
lent potentials, which become overscreened for momenta 
less than a critical value, in agreement with other ana- 
lytical and numerical calculations. 

In conclusion, we believe that our approach provides 
a realistic qualitative description of the paramagnetic 
phase of the diluted magnetic semiconductors. Caution, 
however, should be exercised in applying our calculation 
in the limit of £ approaching unity, where the paramag- 
netic model breaks down and our approximation might 
lead to singularities. 
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APPENDIX A: CORRELATION FUNCTION AT 
THE ORIGIN 

A realistic estimate of the two-particle correlation 
function at the origin in a two-dimensional electron gas 
can be obtained by following an idea developed in Ref. 
|34| . Since the screened Coulomb interaction between elec- 
trons has radial symmetry, a pair of electrons with oppo- 
site spin forms a singlet with spatial wave function that 
depends only on the magnitude of the relative distance 
between them, ^(pi,p 2 ) — ^(p = \pi — fa])- By using 
cylindrical coordinates and dropping the dependence on 
the angle and the perpendicular coordinate we find that 
^f(p) verifies the following Schrodinger equation: 



ti 2 /d 2 * ld#\ 
2m* V dp 2 p dp J 



v(pMp) = m{p) (Ai) 



where V(p) is the effective potential, and m* = m/2 is 
the reduced mass. It is also convenient to define k such 
that E = h 2 k 2 /(2m*). 

We approximate the screened Coulomb potential by 
the potential of an electron that is surrounded by a circle 
of radius r s uniformly filled with screening charge density 
ne = e/(irr 2 ). Outside the circle the charge is zero and 
so the effective potential, 

^2 



V(p) 



P 



±e(£-) + *- 

7T V r s / 7T 

V(p) 



1 



0, 



P^r s 
P^r s 



(A2) 



where E(x] is the complete elliptic integral of the second 
kind. For convenience, we introduce dimensionless vari- 
ables, x = p/(r s aB) and q = kr s aB, where = h 2 /me 2 
is the Bohr radius. The Schrodinger equation becomes: 

d 2 * Id* , T , . 

-TT + —r- + q 2 ^{x) = ; i > 1 

dx z x dx 



d 2 * 

dx 2 



I d* 

x dx 



rl 

X 



-E(x) 



x < 1 



4 

7T 



(A3) 

The general solution for x I is given by ^j(x) — 
AJ (qx) + BN (qx), where A and B are constants, Jo is 
the Bessel function of order zero and Nq the correspond- 
ing Neumann's function. To find the solution inside the 

oo 

circle we make a Taylor expansion of ^(x) = a n x n . 

71=0 

Since we are interested in the ground state of the system, 
its energy is small and so is q. Then, we drop the q 2 term 
from the differential equation, and arrive to the following 
recurrent relation between the a n coefficients: 



n 2 av 



+2 



OO 

E 

rn—l 



(2m- I) 



2" l m! 



2m — 



(A4) 



As a consequence of this recurrent relation, every a n is 
proportional to ao and a f unc tion of r s : a n = &oF n (r s ). 

In order to solve Eq. (A3) we match \&(x) and its 
derivative at x = I : 

*(a; = I) = a Q G(r s ) = AJ (q) + BN (q) 

= 1) = a F(r s ) = AJ Q (q) + BN' (q) (A5) 

OO OO 

where G(r s ) = F n (r s ) and F(r s ) = nF n (r s ). By 

making an expansion in q and keeping the higher orders, 
we arrive to the following relation: 

a G(r s ) = A+ -B[\n(q/2) + C] + 0(q 2 lnq) 

IT 





T 






2 -B 






7T 


Q 







2 \2 



0{q 2 \nq) (A6) 
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TABLE II: Two-particle correlation function at the origin for 
various values of r„. First column displays the values ob- 



tained using Eq (A9), second one the numerical calculation 
by Nagano et al. and last column the interpolation results of 
Polini et al. 



r s 


Eq. ( 


A9 


) 


Nagano et al. 


Polini et al. 




0.5 
1 

2 
5 


0.5 
0.299 
0.199 
0.106 
0.032 


0.50 
0.31 
0.21 
0.13 


0.5 
0.293 
0.204 
0.123 
0.050 



where C is the Euler constant. Therefore, in the limit of 
small momentum 



1 



CtQ oc 



qlnq 



(A7) 



G(r a ) ' G(r.)a'- 
Using the recurrent relation ( |A4j) , we can obtain G(r s ): 

1 1 



G(r s ) = l + r s {- + - 

14 7T 



+ 



OO 

E 



(2m- 1)! 



1 



(2m- l)(2m + 2) 2 
f 0(r 2 s ) ~ (1 + 0.5857r s ) (A8) 



And, since the pair-pair correlation g{p = 0) is propor- 
tional to the square of the wave function ^(p = 0), 



g{p = o) 



2 [1 + 0.5857r s ] 2 



(A9) 



This formula agrees quite well with previous calcu- 
lations, as it can be seen in Table || where the values 
of th e co rrelation function at the origin, obtained using 
Eq. (|A9|), are_cpmpared to the numerical calculation of 
Nagaa© et al.E3 and the most recent estimate by Polini 
et alB 



1 For a review, see G. Prinz and K. Hathaway, Physics To- 
day, 48, (4), 24. 

2 S. A. Crooker, D. A. Tulchinsky, J. Levy, D. D. 
Awschalom, R. Garcia and N. Samarth, Phys. Rev. Lett. 
75, 505 (1995). 

3 H. Munekata, H. Ohno, S. von Molnar, A. Segmiiller, L. 
L. Chang and L. Esaki, Phys. Rev. Lett. 63, 1849 (1989). 

4 H. Ohno, H. Munekata, T. Penney, S. von Molnar and L. 
L. Chang, Phys. Rev. Lett. 68, 2664 (1992). 

5 G. D. Mahan in Many-Particle Physics, chapter 5 (Plenum 
Press, New York, 1990). 

6 F. Pederiva, E. Lipparini and K. Takayanagi, Europhysics 
Lett., 40, 607 (1997). 

7 D. C. Marinescu and I. Tifrea, Phys. Rev. B 65, 113201 
(2002). 

8 G. Niklasson, Phys. Rev. B 10, 3052 (1974). 

9 S. Moroni, D. M. Ceperley and G. Senatore, Phys. Rev. 
Lett. 75, 689 (1995). 

10 D. C. Marinescu and J. J. Quinn, Phys. Rev. B 56, 1114 
(1997). 

11 M. Polini and M. P. Tosi, Phys. Rev. B 63, 045118 (2001). 

12 C. A. Kukkonen and A.W. Overhauser, Phys. Rev. B 20, 
550 (1979). 

13 X. Zhu and A. W. Overhauser, Phys. Rev. B 33, 925 
(1986). 

14 K. S. Yi and J. J Quinn, Phys. Rev. B 54, 13398 (1996). 

15 J. Hubbard, Proc. R. Soc. London, Ser. A 240, 539 (1957); 
243, 336 (1957). 

16 G. Vignale and K.S. Singwi, Phys. Rev. B 31, 2729 (1985) 
and G. Vignale and K.S. Singwi, Phys. Rev. B 32, 2156 
(1985). 

17 We ignore spin-flip processes since the spin-diffusion time is 
several orders of magnitude bigger than typical Coulomb 
scattering times; see I. D'Amico and G. Vignale, Phys. 
Rev. B 62, 4853 (2000). Also the tiny dipole-dipole inter- 



action is ignored. 

18 The correlation between same spin electrons is consider- 
ably smaller than all the other interactions, since same 
spin electrons are kept apart in agreement with the Pauli 
principle. 

19 F. Stern, Phys. Rev. Lett. 18, 546 (1967). 

20 For a point-like charge the density becomes p(q) = 1. 

21 a* B = ao = 0.529 ■ 10 -8 cm for a regular three dimensional 
metal and a* B = 195 ao = 1.04 ■ 10~ 6 cm for GaAs. 

22 J. Moreno and D.C. Marinescu, to be published. 

Note that Vrpa = -, r is equivalent to the 

CRPA{q, oj = 0) 

RPA effective interaction potential, Vrpa(<\,uj = 0), ob- 
tained from equation (Q) by neglecting the local field cor- 
rections. 

24 C. Bulutay and M. Tomak, Phys. Rev. B 53, 7317 (1996). 

25 A. Holas, P.K. Aravind and K.S. Singwi, Phys. Rev B 20, 
4912 (1979); 25, 561 (1981). 

26 O. V. Dolgov, D. A. Kirzhnits and E. G.Maksimov, Rev. 
of Mod. Phys. 53, 81 (1981) . 

27 B. Tanatar and D. M. Ceperley, Phys. Rev. Phys. B 39, 
5005 (1989). 

28 Although there have been attempts to extract the suscep- 
tibility from Monte Carlo simulations, large error bars are 
generally associated with the results. 

29 For the unpolarized systems becomes negative for r s > 
17. 

30 K. S. Singwi, A. Sjolander, M. P. Tosi and R.H. Land, 
Phys. Rev B 1, 1044 (1970). 

31 P. Vashishta and K. S. Singwi, Phys. Rev B 6, 875 (1972). 

32 Y. R. Wang, M. Ashraf and A. W. Overhauser, Phys. Rev. 
B, 30, 5580 (1984). 

33 For example, G^(q — > oo,( — —1) < when r s > 1, but 
G^(q — > oo, C, = —2/3) < when r s — > oo. 



10 



A. W. Overhauser, Can. J. Phys. 73, 683 (1995). 

S. Nagano, K. S. Singwi, S. Ohnishi, Phys. Rev. B, 29, 

1209 (1984). 



M. Polini, G. Sica, B. Davoudi and M. P. Tosi, J. Phys.: 
Condens. Matter 13, 3591 (2001). 



